Frustrated Expectations: Defect Networks and Dark Energy 



O 
O 

(N 

X> 
IX, 

(N 



> ■ 

o : 

in ■ 
(N ■ 

o : 

O ■ 

^ : 

Oh; 

6 ■ 
a: 



13 



P.P. Avelino, 1 ' 2 '0C.J.A.P. Martins, 1 ' 3 'Q J. Menezes, 1 ^ R. Menezes, 1 - 4 '! and J.C.R.E. Oliveira 1 ' 2 ^ 

1 Centro de Fisica do Porto, Rua do Campo Alegre 687, 4169-007 Porto, Portugal 
2 Departamento de Fisica da Faculdade de Ciencias da Universidade do Porto, 
Rua do Campo Alegre 687, 4169-007 Porto, Portugal 
3 Department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences, 
University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, United Kingdom 
* Departamento de Fisica, Universidade Federal da Paraiba, 
Caixa Postal 5008, 58051-970 Joao Pessoa, Paraiba, Brazil 
(Dated: 24 February 2006) 

We discuss necessary conditions for a network of cosmic domain walls to have a chance of providing 
the dark energy that might explain the recent acceleration of the universe. We derive a strong 
bound on the curvature of the walls, which shows that viable candidate networks must be fine- 
tuned and non-standard. We also discuss various requirements that any stable lattice of frustrated 
walls must obey. We conjecture that, even though one can build (by hand) lattices that would 
be stable, no such lattices will ever come out of realistic domain wall forming cosmological phase 
transitions. We provide some simple numerical simulations that illustrate our results and correct 
some misconceptions in the published literature, but a detailed numerical analysis is left for a 
companion paper. 



I. INTRODUCTION 

A number of observational datasets seem to indicate 
that the (local) universe became recently dominated by 
a dark energy component whose gravitational behaviour 
is very similar to that of a cosmological constant. The 
simplest explanation for this would be that a cosmologi- 
cal constant is indeed responsible, but unfortunately the 
corresponding value of the energy density is in violent dis- 
agreement with particle physics expectations, and there- 
fore a considerable effort has been put into finding alter- 
natives. 

A remarkable possibility is provided by topological de- 
fects Q . If our current understanding of particle physics 
and unification scenarios is correct, defect networks must 
necessarily have formed at phase transitions in the early 
universe [3] • Whether or not they are a viable candidate 
to explain the observed acceleration will depend on their 
detailed dynamics. 

It has been shown 3] that a key ingredient is what 
is commonly called frustration: in order to accelerate 
the universe, which necessarily requires a negative pres- 
sure (more specifically w = p/p < —1/3) the networks 
must be frozen in comoving coordinates. In other words, 
they must simply be conformally stretched by expansion, 
and have arbitrarily small velocities. In this limit cos- 
mic strings would have w = —1/3 (which would nearly 
qualify), though it would be hard to reconcile with ob- 
servations. In any case, we will show in Sect. [n]that no 
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string network will ever reach this limit. On the other 
hand, domain walls in the same limit are expected to 
have w = —2/3, which makes them a much more promis- 
ing candidate. 

It goes without saying that domain walls are cosmo- 
logically quite dangerous, and therefore there are very 
stringent bounds on them. These have originally been 
discussed by Zel'dovich et al. Q and later extended in 
various ways 0,0, Q|- However, it must be pointed out 
that these bounds have not been derived from detailed 
studies of domain wall dynamics, but rather are based on 
simple estimates of what this dynamics should be. The 
most glaring example is that it is often implicitly or ex- 
plicitly assumed that there is about one domain wall per 
Hubble volume, which as we shall see when discussing 
the wall dynamics in Sects. IIIII and llVl need not be the 
case — in fact it can't be the case for a frustrated network. 
There are also more detailed but purel y ph enomenologi- 
cal analyses of this type of models [aliJIlfl- (The second 
of these in fact claims that under certain fine-tuned con- 
ditions on cosmological parameters frustration may not 
be needed.) At a qualitative level these provide useful 
indications of allowed ranges for cosmological and model 
parameters, but again the analyses are somewhat mis- 
led by incorrect assumptions (made by others) on the 
dynamical properties of the networks. 

By the same token, there is also no through study of 
the conditions under which domain wall networks may 
become frustrated. Some relatively low resolution 2D 
numerical simulations exist [Til IT^ | , with mixed results. 
On one hand |Tlj finds that even in non-trivial models 
(Zn models with moderate values of N) annihilation pro- 
cesses can be fairly efficient, so that the networks show 
no tendency towards frustration. On the other hand, 
for fairly large N and suitable choices of initial condi- 
tions yj| finds some hints of the possible formation of a 
hexagonal-type lattice, though given the small dynami- 
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cal range of the simulations the results should be consid- 
ered inconclusive. The fact that the wall velocities are 
high (assuming they are being correctly measured [lflj ) 
also raises some doubts about the interpretation of these 
results. In another context, a study dealing only with 
compact dimensions |l4| also finds no frustration. 

Last but not least, there have been recent attempts to 
construct (by hand) plausible domain wall lattices, and 
then study their stability 0, 0, . Again these can 
provide interesting and suggestive results, but of course 
they fail to address the key issue of if and how such lat- 
tices may emerge from realistic initial conditions for do- 
main wall forming phase transitions. 

The current report is the first in a series that aims to 
put the study of domain wall frustration (or lack thereof) 
on a firmer basis. We shall start by making use of a re- 
cently developed analytic model for domain wall evolu- 
tion (the basics of which are explained in Sect. to 
derive, in Sect. IIVI very general necessary conditions for 
frustration. The apparently innocuous requirements of 
domination around the present epoch and small velocities 
are enough to show that any candidate networks must 
be quite un-naturally fine-tuned (in a precisely quantifi- 
able way). After this 'top-down' result we'll proceed in 
Sect. [V]to discuss an equally important 'bottom-up' one, 
namely some requirements that any stable lattices must 
obey. A related analysis is made in [l^, where energetic 
arguments are discussed for a specific model. The discus- 
sion will concentrate on the 2D case (which is physically 
much simpler than 3D) but some results are expected to 
be generic. Finally, in Sect. IVII we present a summary 
of our results. We will also provide a few very simple 
numerical illustrations in Sect. [V] though we leave an 
extensive numerical study for a companion paper. 



II. COSMIC STRINGS— A NON-STARTER 

Before we focus on domain wall networks, it is worth 
starting with a brief note of the case of cosmic strings, 
which are sometimes claimed to be able to accelerate the 
universe if they are the dominant energy density com- 
ponent. Recall that the equation of state of a cosmic 
string network depends on its root-mean squared (RMS) 
velocity as 



(1) 



so in the limit v s — its correlation length would behave 
as 



L oc a oc t . 



(2) 



Notice that there is a crucial point here which is often 
overlooked: because one knows that there is a sizable 
matter component S7 m0 ~ 0.3 whose equation of state 
is w m = 0, then in order to accelerate the universe one 
would need the dark energy component to have Wd < 



—0.5. Be that as it may, it has been shown ,201 that when 
a string network dominates the universe, the behaviour 
of its correlation length and RMS velocity is 



L — Qt cx a 
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where /i is the string mass per unit length, and 

v = const , 
while the scale factor evolves as 

a cx t 1 ' 1 ^ 2 . 



(3) 



(4) 



(5) 



(6) 



As expected, the naive solution (0 is recovered for v = 0. 
Note that the string correlation length grows as fast as 
allowed by causality (this is always the attractor for the 
evolution of a string network), but the network is not 
conformally stretched (which would be the case only if 
v = 0). Since the string velocity does not vanish, the 
correlation length grows faster than the scale factor. On 
the other hand the effect of the non-zero velocities is to 
make a string-dominated universe expand more slowly 
than one might naively expect, since some of the energy 
of the string network is lost to velocity redshift. 
Numerically it is found [2(j that 



0.17 



(7) 



for non-intercommuting strings. The velocity is expected 
to be smaller, but still non-zero, for the case of entangled 
string networks (unlike in the former case, it is not easy 
to numerically determine this value). Ancillary evidence 
for these results has also been recently discussed in [2l| . 
Hence in a string-dominated universe the string velocity 
never becomes arbitrarily small, and a string- dominated 
universe will not in any circumstances frustrate or accel- 
erate the universe. 

Physically, this result is to be expected: string evolu- 
tion is a non-equilibrium, irreversible process (cf. string 
evolution in a contracting universe 0, H), and their 
dynamics naturally leads them towards relativistic ve- 
locities. In other words, non-trivial mechanisms must 
be active if they are to remain non-relativistic, let alone 
freeze completely. Note that friction due to particle scat- 
tering, which is the simplest such mechanism, can only 
have a limited effect. To some extent these points will 
also relevant for the case of domain wall networks. 



III. DOMAIN WALL EVOLUTION 

Domain walls arise in models with spontaneously bro- 
ken discrete symmetries 0, . The simplest example is 
that of a scalar field cf> with the Lagrangian 



^ a - V{<j>) 



(8) 
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where the potential V(4>) has a discrete set of degenerate 
minima, say for example 



(9) 



V(t)=V ^ Q -l 
By varying the action 

S = I dt { d 3 Xy^jC 



(10) 

with respect to <f> we obtain the field equation of motion 

dV 



d 2 t dcj> 



W = — 



(11) 



where V is the Laplacian in physical coordinates and 
H = (da/dt)/a is the Hubble parameter. 

In many cosmological contexts of interest, one can 
neglect the domain wall thickness when compared to 
its other dimensions, and thus treat the wall as an in- 
finitely thin surface. With this assumption, its space- 
time history can be represented by a 3D world-sheet 
x t* = x fifsa^ a = o, 1, 2. A new action can then be easily 
derived |l|. l n the vicinity of the world-sheet a conve- 
nient coordinate choice is the normal distance from the 
surface. Noticing that in the thin wall limit all fields in 
the Lagrangian should depend only on this normal coor- 
dinate, and integrating out this dependence, one finds 



S = -a d\^j, 



where 



lab — 9 [ivX , a x ,b 



(12) 



(13) 



is the world-sheet metric, with the obvious definition 7 = 
det("f ao ), and a is the mass per unit area of the wall. 
Notice that this action is proportional to the 3- volume of 
the wall's world-sheet, and hence is clearly the analogue 
of the Goto-Nambu action for strings. 

In analogy with the velocity-dependent one-scale 
model for cosmic strings 0, |2^, |2(j, one can obtain a 
one-scale model for domain wall evolution. This was de- 
rived and tested against numerical simulations in |13| so 
here we simple state the results we shall be using. Let us 
define a characteristic length scale, 



L= a - 
P 



(14) 



which is directly related to the average distance between 
adjacent walls measured in the frame comoving with the 
expansion of the universe. One can then show that its 
evolution equation is as follows 



^ = (l + 3v 2 )HL + c w v. 



(15) 



where c w is a phenomenological parameter measuring the 
efficiency of energy losses from the wall network. Here v 



is the RMS velocity of the walls, which in turn evolves 
according to 



* (1 
dt y 



h 
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— 3Hv 



(16) 



with k w being another phenomenological parameter re- 
lated to the typical curvature of the walls — see |2(j for a 
thorough discussion of the analogous parameter for cos- 
mic strings. These therefore provide a phenomenological 
model for domain wall evolution with two free parame- 
ters, that one can measure from high-resolution numeri- 
cal simulations 0, |2jJ, |2£| . 

It is easy to see that, just as for cosmic string networks, 
the attractor solution to these evolution equations corre- 
sponds to a linear scaling solution 



L = et. 



v = const . 



(17) 



Assuming that the scale factor behaves as a cx t a the 
detailed form of the above linear scaling constants is 



3a(l — a) 



1 — a 
3a k u 



(18) 



(19) 



As in the case of cosmic strings [2(j , an energy loss mech- 
anism (that is, a non-zero c w ) may not be needed in order 
to have linear scaling. Note that this means that having 
non-intercommuting domain walls is by no means suffi- 
cient to ensure a frustrated wall network! Note, however, 
that the linear scaling solutions are physically very dif- 
ferent for cosmic strings and domain walls. In the case 
of cosmic strings, in the linear scaling phase the string 
density is a constant fraction of the background density, 
whereas in the case of domain walls we have 



Pw , 

— cx t . 
Pb 



(20) 



so the wall density grows relative to the background den- 
sity, and will eventually become dominant. 

There is in general, however, an effect which we have 
neglected thus far. At early times, in addition to the 
damping caused by the Hubble expansion, there is a fur- 
ther damping term coming from friction due to particle 
scattering off the domain walls. Phenomenologically, it 
can be shown [lj that its effect can be adequately de- 
scribed by a frictional force per unit area 



-7V . 



(21) 



In the above we have defined a friction length scale which 
we write as 



(22) 



/3T 4 
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where T is the photon temperature and (3 is a parame- 
ter counting the number of species interacting with the 
domain walls. Note that it is clear that they can't inter- 
act strongly with the photons, baryons and dark matter, 
otherwise they would be ruled out due to the strong cos- 
mological signatures left over on the cosmic microwave 
background and/or large scale structure El |M E|. 

Just like in the case of cosmic strings |24l I2BI l26| 
one can modify the evolution equations of our one-scale 
model to account for this extra friction term. They be- 
come 

—v 2 + c w v (23) 



IV. 



WALL NETWORK PROPERTIES 



dt 



where we have defined a damping length scale 

h= 3H+ r f m> 

which includes both the effects of Hubble damping and 
particle scattering. It is instructive to compare the im- 
portance of both of these terms. At early times particle 
scattering is always dominant (except if it is completely 
absent, that is if the network is totally non-interacting), 
but the friction length scale grows faster that the Hubble 
length and so eventually particle scattering will become 
sub-dominant. The epoch at which the transition occurs 
depends both on the parameter (3 and on the mass scale 
of the walls, which we can characterise by the scale 77 of 
symmetry breaking phase transition that produced them, 

cr ~ if . (26) 

One can easily find, for j3 ~ 1, that friction domination 
will end at the epoch of equal matter and radiation den- 
sities for a symmetry breaking scale 



Veq ~ {mpiTj q ) 



1/3 



3GeV 



and will end around the present day for 



T, 



<:'<! 



1/6 



1 MeV , 



(27) 



(28) 



where To and T eq are respectively the CMB temperatures 
at the present day and at the epoch of equal matter and 
radiation densities. Note that the dependence on the pa- 
rameter (3 is weak, and in any case the parameter is quite 
tightly constrained (it can't be much larger than unity). 
This coincides with the Zel'dovich bound Q, though we 
emphasise that the derivation is different. In particu- 
lar, the classical derivation implicitly assumes the linear 
scaling solution, which as we saw above is not generically 
the case: one defect per correlation volume does not nec- 
essarily imply one defect per Hubble volume (this point 
has also been made in [la]). In the following section we 
shall see that tight constraints apply to the parameters 
of any domain wall network if we assume that it becomes 
frustrated and is starting to accelerate the universe. 



If a domain wall network is to provide the dark energy 
suggested by observations, there are some obvious and 
unavoidable requirements. Firstly, such a network must 
be dominating the energy density of the universe around 
the present day, so its energy density must be of the order 
of the critical density, 



Pu 



1 



(29) 



which provides us with a unique relation between the 
energy scale of the defects and the present correlation 
length, namely 



if 



T T eq 



(30) 



But there is a further constraint on the physical cor- 
relation length today. The dark energy should be ap- 
proximately homogeneous and isotropic on cosmological 
scales or otherwise that would result in strong (unob- 
served) signatures on the cosmic microwave background. 
So the product 



LqHo 



f_jL_y 

\30MeVJ 



(31) 



,30MeV; 

must be much smaller than unity. If we say we need 

L < IMpc << H- 1 , (32) 



then we find again 



77 < IMeV . 



(33) 



Now, recall that the averaged equation of state of a 
domain wall network is given by jlflj 



w v 



2) 



(34) 



where now v w is the averaged RMS velocity of the do- 
main wall network. So in order to accelerate the universe 
with an equation of state in agreement with observations 
the wall velocities must necessarily be quite small. In 
general, this means that the network should be friction- 
dominated, which takes us back to the bound l|28|) , which 
coincided with (|33ll . However, during friction domina- 
tion, one also has |13| 



Lv ~ k£f ~ fc^j 



(35) 



from which we get at the present day, using the relation 
(EOJ), 



T 

T ' 



(36) 



and since the velocity must be much less than unity, 

fc«l(T 4 . (37) 
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The only way to avoid this bound is to enforce that 
one has a model where there is no friction due to particle 
scattering (effectively if = oo). In such a case iPDESjl 
still apply, but instead of being friction-dominated all the 
way through to the present day the network would be in 
the linear scaling regime. In that case, (|35|) is replaced 
by 



V. 



WALL LATTICE PROPERTIES 



Lv ~ kt ~ — , 
Ho 



(38) 



and again substituting Lq and requiring that the velocity 
is much smaller than unity we have 

k < L H « 1(T 4 . (39) 



Hence we see that the curvature of the domain walls must 
unavoidably be very small. Note that for the case of or- 
dinary cosmic strings |25l l2fi| | k is a parameter depending 
on the string velocity, whose value increases and closely 
approaches unity in the limit of small velocities, and it is 
expected that the same happens for the simplest domain 
wall models. So the only possibly realistic candidates 
are non-standard networks, that is those with junctions 
where N > 2 walls intersect. 

It is also important to emphasise that for any realistic 
network that is becoming the dominant energy compo- 
nent of the universe around the present day, assuming 
a constant equation of state is generically not a good 
approximation. Because the network's equation of state 
depends on the network velocity (see eqn. 1341) , it will only 
be constant if the network is in a linear scaling regime (in 
which case w > —2/3) or in the asymptotic limit v = 0. 
In all other cases w should be time-varying. Assuming 
w = const may therefore be a poor approximation. This 
is one of the reasons why currently existing phenomeno- 
logical analyses of this class of models [a H3 are omv 
at best qualitatively accurate, and a more realistic study 
is called for. 

We note that the velocity dependent one scale model 
for domain walls does not fully take into account the fea- 
tures of more complex models, since in that case there 
might be added contributions (in particular those related 
with the dynamics of the junctions). However, one can 
argue that the dynamics of these junctions will be re- 
flected in the evolution of curvature parameter which in 
our model is simply a phenomenological parameter char- 
acterising domain wall curvature. Hence, we do expect 
that our simple one scale model model to also be valid (at 
least qualitatively) in the context of these more complex 
configurations. Be that as it may, in the next section we 
shall discuss some properties of such wall network lat- 
tices. 

Finally, we note that one of the assumptions in the 
model is that the domain wall mass per unit area, tr, is 
fixed. It is certainly possible to envisage more complex 
models with different types of walls with various masses 
per unit area. However, we do not expect our conclusions 
to be modified even in these more complex scenarios. 



In the previous section we have shown that, in order to 
be able to provide the dark energy, the curvature of the 
domain walls would have to become very small by the 
present time. In this section we shall look at geometrical 
properties of polyhedrons in order to investigate if equi- 
librium flat domain wall configurations can be the natural 
result of domain wall network evolution in 2-dimcnsions. 

The Poincare formula relating the number of polyhe- 
dron vertexes (V), faces (F), and edges (E) of genus g 
surfaces has the form 



V-E + F = 2-2g. 



(40) 



We shall assume periodic boundary conditions on a two 
dimensional square box and consequently one is effec- 
tively considering a surface of genus equal to unity. Hence 
eqn. f4"U|) becomes 



V- E+F = . 



(41) 



We note that our assumption of periodic boundary con- 
dition does not affect our conclusions as long as the size 
of the box is big enough. 

Let us start by considering the case in which the num- 
ber of edges of each polygon, x, and the number of edges, 
d, meeting at a vertex are fixed. Let us denote the num- 
ber of polyhedron faces by F — N x . The number of 
polyhedron vertices is V = N x x/d since each polygon has 
N x vertexes but each one of them is shared with d — 1 
other polygons. Also the number of polyhedron edges is 
equal to E = N x x/2 since each polygon has N x edges 
but each one of them is shared with another polygon. 
Consequently, in this case eqn. l|4*T|l becomes 



N x 1 + - - 



= 



(42) 



This equation has the following solutions (x = 6, d = 3), 
(x = 4, d = 4), (x = 3,d = 6). These are the well 
known hexagonal type lattices with odd Y-type junc- 
tions, square lattices with even X-type junctions and tri- 
angular lattices with even V-type junctions in 2 dimen- 
sions. 

However, in general we do not expect that all the poly- 
gons have the same number of edges. Causality con- 
straints mean that such configurations could not arise 
directly out of a cosmological phase transition, though 
they could conceivably be generated dynamically (more 
on this below). Hence, let us consider the more interest- 
ing situation in which d is fixed but x is not. In this case, 
it is straightforward to show that 



2d 



EZiN x d-2 



(43) 



It is also easy to show that (x) = 6 if d = 3, (x) =4 
if d = 4, (x) = 3 if d = 6 and (s) ->■ 2 if d -> oo. We 
will show, in the following discussion, that these simple 
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FIG. 1: An illustration of the decay of an unstable X-type 
junction into stable Y-type junctions. The decay is energeti- 
cally favourable since it leads to a reduction of the total length 
of the walls (if all have the same tension). 



geometrical considerations will be relevant for evaluating 
the potential of domain walls as a dark energy candidate. 

There are no polygons with two edges. However, the 
domains in a realistic domain wall network will not in 
general have straight edges and consequently two edge 
domains are possible. However, these domains will be 
unstable and collapse due to the domain wall curvature 
independently of the number of elements meeting at each 
junction. 

Here, we are implicitly assuming that the energy asso- 
ciated with the junctions is negligible which in practise 
means that they are free to move. This is a reasonable 
assumption, at least for the purposes of the present dis- 
cussion. If it were not the case, one would have to take 
into account the contribution of the junctions when cal- 
culating the equation of state associated with the domain 
wall network, and consequently w = p/ p would necessar- 
ily be greater than —2/3 even for a fully static configu- 
ration. Such networks would hardly be compatible with 
observational bounds. 

It is also straightforward to show from a local stability 
analysis that three, four and five edge domains will also 
be unstable if only Y-type junctions {d = 3) occur in 
a given model. This is a particularly interesting case 
since one can show using local energy considerations that 
in models with more than two vacua, if all the domain 
walls connecting the various vacua have equal energies 
then only Y-type stable junctions would form (see Fig. 

Higher-order junctions are unstable and very quickly 
decay into Y junctions. A good illustration of this point 
is the 'pentahedral' model discussed in [l^. The author 
erroneously claims that this a candidate for frustration 
with X-type junctions, when in fact it will form Y-type 
junctions. This can easily be checked numerically, and 
so can the fact that even if one constructs (by hand) a 
box with X-type junctions these will quickly decay into 
Y-type junctions. 




FIG. 2: An illustration of the collapse of three (top) and four 
(down) edge domains with Y-type junctions. The collapse is 
energetically favourable since it leads to a reduction of the 
total length of the walls. 




FIG. 3: An illustration two different six-edged polygons with 
Y-type junctions of walls with same tension. Both configura- 
tions have the same energy. 

Fig. [2] shows various polygons formed by walls with 
equal tension. Note that in both cases, x = 3 (top) and 
x — 4 (bottom), the total length of walls decreases. Con- 
sequently, the polygons will tend to collapse thus min- 
imising their potential energy. On the other hand in Fig. 
[3] for x = 6, the length remains constant and both config- 
urations have the same energy. Hence, given that d = 3 
implies (x) = 6, the only possible equilibrium configu- 
ration with only Y-type junctions is a hexagonal type 
lattice (otherwise unstable two, three, four and five edge 
domains would occur). 

Although hexagonal lattices with Y-type junctions al- 




FIG. 4: An illustration of the collapse of a four-edged polygon 
in the case where two of the surrounding domains are on the 
same vacuum state. The collapse leads to a reduction of the 
number of edges of contiguous domains. 
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FIG. 5: An illustration of the collapse of a four-edged polygon 
in the case where all the surrounding domains are on different 
vacuum states and the domain walls all have the same tension. 
Again, the collapse leads to a reduction of the number of edges 
of contiguous domains. 



low for locally confined energy conserving deformations, 
Hubble damping may prevent the collapse of such config- 
urations (even if they are perturbed) and consequently 
one should not completely discard them on this basis. 
However, we do not expect that hexagonal lattices will 
be an attractor for the evolution of domain wall network 
simulations. Although a four edge domain is unstable 
in this context, its collapse will modify the properties of 
its four contiguous domains with Ej, Ejj, Em and Eiy 
edges. If we assume that the first two will join, as hap- 
pens in the example given in Fig0] the resulting domain 
will have Ej + Ejj — 4 edges while the other two will have 
Em — 2 and Ejy — 2 edges after the collapse of the four 
edge domain. The production of three hexagons as a re- 
sult of the collapse of a four edge domain is improbable 
since it would require Ej+Ejj = 10 and Em = Ejy = 8. 

Of course, in a model with a very large number of vacua 
the probability that two nearby triple wall junctions will 
annihilate can be made arbitrarily small. In fact as the 
number of vacua increases, the probability that the con- 
figuration shown in Fig. [3] occurs becomes much greater 
than the one given in Fig. 0] Still, the collapse of un- 
stable domains with two, three, four and five edges will 
always result in a decrease in the number of edges of some 
of the contiguous domains. Again, we do not expect that 
a domain wall network in two dimensions will naturally 
evolve towards a hexagonal lattice from realistic initial 
conditions. 

Also, no equilibrium configurations exists with d > 6. 
This means that if we started with a domain wall net- 
work with d > 6 unstable two edge domains would neces- 
sarily be present and consequently the number elements 
meeting at a junction will often have to be reduced to a 
number smaller than 6. Hence, we anticipate that odd Y - 
type and/or even X-type junctions will be generic even 
in models where the number of elements meeting at a 
junction is allowed to be greater than four. 

We have seen that in a model where all the domain 
walls connecting the various vacua have the same tension, 
it is not expected that a frustrated domain wall network 
will naturally occur. However if we relax this assump- 
tion we are adding a different source of instability since 
the walls with higher tension will tend to collapse thus 
increasing the dimensionality of the junctions which, in 




FIG. 6: An illustration of the collapse of two Y-type junctions 
into one X-type junction. The thickness of the traces indicates 
the tension strength. If the thick wall has a tension larger 
than twice that of the lower tension ones the collapse must 
happen. 




FIG. 7: The evolution of a perturbed square lattice with even 
X-type junctions in a matter-dominated universe. The top left 
panel is the initial configuration. From left to right and top to 
bottom panels the horizon is approximately 1/256, 1/20, 1/10 
and 1/5 of the box size respectively. The lattice stabilises in 
the right bottom panel configuration. 

turn, will lead to the production of further unstable two 
edge domains. This is clearly illustrated in Fig. which 
shows the collapse of two Y-type junctions into one X- 
type junction which must necessarily occur if the thick 
wall has a tension larger than twice that of the lower 
tension ones. 

Hence, although we have not presented a rigorous 
formal proof, we conjecture that it is unlikely that 2- 
dimensional domain wall network evolution from realistic 
domain wall forming phase transitions will ever produce 
a frustrated network. Our analysis can only be fully ap- 
plied to domain walls networks in 2-dimensions. How- 
ever, some of our results can also be generalised to 3- 
dimensions, at a cost of a much greater complexity. It is 
not clear that the increase in the number of degrees of 
freedom which occurs when we go from 2 to 3-dimensions 
will help frustration and consequently we believe that our 
simpler analysis in 2-dimensions is a crucial step in ac- 
cessing domain walls networks as dark energy candidates. 

Finally, we illustrate our results with two simple nu- 
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FIG. 8: The evolution of a perturbed hexagonal lattice with 
odd Y-type junctions in the matter-dominated epoch. The 
top left panel is the initial configuration. From left to right 
and top to bottom panels the horizon is approximately 1/256, 
1/10, 1/5 and 2/5 of the box size respectively. The lattice 
stabilises in the right bottom panel configuration. 

merical examples. The stability of simple examples of 
triangular, hexagonal and square domain wall lattices in 
2D was studied by analytic means in 

EE EL specif- 
ically by looking at their macroscopic elastic properties. 
The authors assess the stability of various lattice config- 
urations and find that some of them are stable. Here, we 
confirm using numerical simulations that there are cer- 
tain lattice configurations which are indeed stable even 
against large deformations. In Fig[7|we show a field the- 
ory numerical simulation in the matter dominated epoch 
for a perturbed square lattice with even X-type junc- 
tions. As expected the network evolves towards the min- 
imum energy equilibrium configuration. The evolution 
of a perturbed hexagonal lattice with odd Y-type junc- 
tions in the matter dominated epoch is shown in FiglHI 
As it was mentioned before this is the only possible equi- 
librium configuration if only Y-type junctions are per- 
mitted. Although in jl5l Ha . the authors claim that 
an hexagonal lattice with Y-type junctions is unstable, 
the expansion of the Universe will damp the domain wall 
velocities which may drive the network towards an equi- 
librium configuration (as shown in Fig|SJ). We have per- 
formed a similar simulation in Minkowski space and in 
that case the configuration is always unstable. This is 
to be expected since no damping mechanism operates in 
this case. These results clearly indicate that the cru- 
cial question is not the existence of specific stable lattice 
configurations but whether any of these can be the natu- 
ral result of domain wall network evolution from realistic 
initial conditions. A more detailed numerical analysis of 
these issues will be left for a forthcoming publication. 
VI. CONCLUSIONS AND OUTLOOK 

We have studied necessary conditions for a network of 
cosmic domain walls to become frustrated, and thereby 



have a chance of accounting for the dark energy. We have 
made use of simple analytic tools to derive a strong bound 
on the wall curvature, which implies that any candidate 
network must be considerably fine-tuned. We have also 
considered various simple lattice properties in two spatial 
dimensions, and used energy considerations to obtain a 
number of requirements that any stable lattice of frus- 
trated walls must obey. Although the extrapolation of 
some of the results of the latter analysis to the case of 
three spatial dimensions is non-trivial, it is clear that 
some general trends will remain. 

Several key points emerge from our analysis, each of 
which will be considered in detail in forthcoming work. 
Firstly, there is some model dependence involved, which 
is to be expected given the whole zoo of existing mod- 
els. Some of these have been studied in the past, but 
others remain explored. Secondly, energy considerations 
are the main driving mechanism for the evolution of wall 
networks with junctions, though topological arguments 
also play a role. An example of the power of these mech- 
anisms is the simple yet crucial result that in models 
with more than two vacua, if all the domain walls con- 
necting the various vacua have equal energies then only 
Y-type stable junctions would form. As we have pointed 
out above, this point has been overlooked in some of the 
existing literature. 

Last but not least, the difference between designer con- 
ditions (that is, lattice configurations built by hand) and 
conditions that might possibly arise as the outcome of 
domain wall forming phase transitions cannot be overem- 
phasised. In this regard, even the presence or absence of 
expansion can have non-trivial effects. Our results lead 
us to conjecture that, even though one can build (by 
hand) lattices that would be stable, no such lattices will 
ever come out of realistic domain wall forming cosmo- 
logical phase transitions. If so, then defect networks are 
ruled out as an explanation for the dark energy. We shall 
explore this conjecture, in particular through detailed nu- 
merical analysis, in a companion paper. 
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